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Abstract — We analyze the effect of feedback delay and channel 
estimation errors in a decode-and-forward (DF) relay selection 
scenario. Amongst all the relays that decode the source informa- 
tion correctly, only one relay with the best relay-to-destination 
(R — > D) channel quality is selected. More specifically, the 
destination terminal first estimates the relay-destination channel 
state information (CSI) and then sends the index of the best 
relay to the relay terminals via a delayed feedback link. Due 
to the time varying nature of the underlying channel model, 
selection is performed based on the old version of the channel 
estimate. In this paper, we investigate the performance of the 
underlying selection scheme in terms of outage probability, 
average symbol error rate (ASER) and average capacity. Through 
the derivation of an ASER expression and asymptotic diversity 
order analysis, we show that the presence of feedback delay 
reduces the asymptotic diversity order to one, while the effect of 
channel estimation errors reduces it to zero. Finally, simulation 
results are presented to corroborate the analytical results. 



Index Terms — Selection cooperation, 
feedback delay, channel estimation error. 
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I. Introduction 

IT has been demonstrated that cooperative diversity pro- 
vides an effective means of improving spectral and power 
efficiency of wireless networks as an alternative to MIMO 
systems (TJ, Q. The main idea behind cooperative diversity 
is that in a wireless environment, the signal transmitted by 
the source (S) is overheard by other nodes, which can be 
defined as "relays" [3]. The source and its partners can then 
jointly process and transmit their information, thereby creating 
a "virtual antenna array", although each of them is equipped 
with only one antenna. It is shown in ||4], that cooperative 
diversity networks can achieve a diversity order equal to the 
number of paths between the source and the destination. 
However, the need of transmitting the symbols in a time 
division multiplexing (TDMA) fashion reduces the maximum 
achievable capacity improvements. Additionally, due to power 
allocation constraints, using multiple relay cooperation is not 
economical. To overcome these problems, relay selection is 
proposed to alleviate the spectral efficiency reduction caused 
by multiple relay schemes and to moderate the power alloca- 
tion constraints IB, GQ, 0- 
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In [6], Bletsas et al. proposed a fast selection algorithm rela- 
tive to the coherence time of the channel. For each relay node, 
they assume timers that are inversely proportional to either 
the harmonic mean or a max — min function of the back-to- 
back channels. They study the performance of their selection 
cooperation scheme in the amplify-and-forward mode (S-AF) 
in terms of the outage probability. The outage probability and 
the average symbol error rate (ASER) expressions of a relay 
selection with decode-and-forward (S-DF) scenario have been 
studied in [8] and [9], respectively. In 1101 . Beres and Adve 
analyze relay selection in a network setting, i.e., multi-source 
networks and introduced a closed-form formula for the outage 
probability of a single source single destination network with 
multiple relays in the DF protocol. In ifTTl . lfl2l the nearest 
neighbor selection for the DF mode is proposed. In their 
work, the best relay is simply the spatially nearest relay to 
the base station. In J5], Q, the relay whose path introduces 
the maximum SNR is selected, in the amplify-and-forward 
(AF) scheme. In J7], the improvement of S-AF over the all 
participate cooperation scheme in the AF mode (AP-AF) is 
discussed in terms of symbol error rate. The ASER of S-AF 
and AP-AF is discussed based on an approximation of the 
cumulative density function of the received SNR. In |5|, Zhao 
and Adve study the outage probability of S-AF scheme and 
analyze the performance improvement due to S-AF. They also 
show dramatic improvements of S-AF over AP-AF in terms of 
capacity. In |[T3l the authors discussed the outage probability, 
ASER, ergodic capacity and outage capacity in S-DF protocol 
and compared their results with selection relaying (SR). In 
fl4l . Ikki and Ahmed derive outage probability, ASER, and 
capacity expressions for S-AF and illustrate the improvement 
of S-AF over AP-AF. It is shown that the relay links in S-AF 
as well as AP-AF introduce the same diversity order of M, 
where M is the number of relays. 

Related work and contributions: The effect of feedback 
delay on the outage probability of antenna selection in a 
MIMO communications system is discussed in fl5l . Relay 
selection in the setups in |6|, II81- 1T21 . assume the perfect 
channel state information (CSI) knowelge. However, in prac- 
tical scenarios, the communication links are not known and 
have to be estimated. Since channel estimation is necessary 
for the selection procedure, channel estimation error due to 
imperfection of the estimator is inevitable. The effect of 
channel estimation error in selection cooperation in the AF 
mode is studied in lfl6l . flTl . In fl6l . we analyze the effect of 
channel estimation error on the outage probability of selection 
cooperation in the AF mode. Capacity of selection cooperation 
with imperfect channel estimation is studied in [17], |18|. The 
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impact of channel estimation error on the ASER performance 
of distributed space time block coded (DSTBC) systems is 
investigated in |fl9l , assuming the amplify-and-forward pro- 
tocol. Building upon a similar set-up, Gedik and Uysal [ 20 1 
extend the work of |[T9l to a system with M relays. In l2D . 
the symbol error analysis is investigated for the same scenario 
as in (20). 

In a practical relay selection scenario, the destination is 
responsible for estimating the CSI and performing relay se- 
lection. The index of the selected relay is then fed back to 
all the relays. Due to the time varying nature of the fading 
channels, which is function of Doppler shift of the moving 
terminals, the CSI corresponding to the selected relay is time 
varying. In this case, relay selection is performed based on 
outdated CSI. Hence, the selection scheme may not yield the 
best relay. 

Although there have been research efforts on conventional 
MIMO antenna selection with feedback delay and channel 
estimation errors (see for example [ 15] and references therein), 
only a few isolated results have been reported in the context 
of cooperative communications. In ||23l , Vicario et al. analyze 
the outage probability and the achievable diversity order of an 
opportunistic relay selection scenario with feedback delay. In 
this paper, we investigate the impact of feedback delay and 
channel estimation errors on the performance of a cooperative 
diversity scheme with relay selection. To the best of our 
knowledge, this is the very first paper that investigates relay 
selection with outdated channel estimates. Assuming imperfect 
channel estimation and feedback delay, our contributions are 
summarized as follows: 

> We derive exact ASER and outage probability expres- 
sions for relay selection with decode-and-forward (S-DF) 
relaying in the presence of feedback delay and channel 
estimation errors. 

> We derive a lower bound on the average capacity in the 
presence of feedback delay and channel estimation errors. 

> We demonstrates that the asymptotic diversity order is 
reduced to one with feedback delay and reduced to zero 
in the presence of channel estimation error. 

• We present a comprehensive Monte Carlo simulation 
study to confirm the analytical observations and give 
insight into system performance. 

The rest of the paper is organized as follows: In Sec. [II] 
we introduce our system setup and the selection strategy. In 
this section the channel estimation error as well as the delayed 
feedback link models are illustrated. In Sec. [Ill] we investigate 
the outage probability of the system and develop a closed-form 
expression for outage probability. In Sec. IIVI we study the 
ASER performance of the system in detail and derive exact 
expression for the ASER. In Sec. [V] we analyze the asymptotic 
order of diversity in S-DF. In particular, we investigate the 
effect of channel estimation errors and feedback delay on 
the diversity of the system. In Sec. [VI] we introduce a lower 
bound for the instantaneous capacity and derive a closed-form 
analytical expression for the average lower bound capacity. 
Simulation results are presented in Sec. IVHI and the paper is 
concluded in Sec. [XI] 
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Fig. 1. System model for decode-and-forward relay selection with delayed 
feedback. 

II. SYSTEM MODEL 

Fig.Q]shows the selection cooperation system model studied 
in this paper. We consider a multi-relay scenario with M 
relays. We assume that the relay R m , m — 1,...,M, the 
source S, and the destination D are equipped with single 
transmit and receive antennas, respectively. In our system 
model, we ignore the direct transmission between the source 
and its destination due to shadowing. h sm and h md represent 
the channel fading gains between S — > R m and R m — > D, 
respectively. Assuming a half duplex constraint, the data trans- 
mission is performed in two time slots. In the first time slot the 
source terminal transmits its data to all potentially available 
M relays. After, receiving the source signal via independent 
channels, all the relays R m , m = 1,2, ...M, decode their 
received signal, and check whether the transmitted signal is 
decoded correctly or not. This can be done via some ideal 
cyclic redundancy codes (CRC) |28|, which are added to the 
transmitted information symbols. We define the decoding set 
T>(s) as the set of relays that decode the transmitted signal 
correctly. Clearly only those relay nodes with a good source 
to relay channel can be in the decoding set T>(s). In the 
second time slot, the best relay that satisfies an index of merit 
participates in the transmission and broadcasts its decoded 
symbol towards the destination. 

In this paper, conditioned on the decoding set T>(s), the 
selection function selects the relay with the best downlink 
channel i.e., 

m* = arg max {j m d} 



arg max {P\h md \ 2 }, 

mGV(s) 



(1) 



where j m d is the received SNR from the m th path. We assume 
that the source has a power constraint of P Joules/symbol 
and similarly each relay node in T>(s) can potentially transmit 
its information with P Joules/symbol, and the receiver noise 
power is 1. 

A. model of channel estimation error 

All channels are assumed to be independent and identically 
distributed (i.i.d), with h ~ CAf(Q, <r?). A block fading scheme 
is considered where the channel realizations are assumed to 
be constant over a block and correlated across blocks. The 
correlation coefficient between the k th and the (k + i) th block 
is p f = J (2irf d Ti), where f d is the Doppler frequency, J (-) 
is the zeroth order modified Bessel function of first kind and 
T is the block duration. Minimum mean square error (MMSE) 
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estimators are employed at the relay nodes. We assume that 
the number of training symbols in each frame is sufficiently 
less than the data symbols, so that increasing the training 
symbols' power would not lead to an increase in the average 
power; hence, we assume that the training symbols' power is 
independently adjustable. Let the true channel be h and the 
estimated channel be h, then h and h are related by [22] 



h = h + e, 
h = p e h + u, 



(2a) 
(2b) 



where e and u are the channel estimation errors, h, e and u, are 
all zero mean Gaussian random variables with variances cr? al 
and a\. p e = (c^/erf ) is the correlation coefficient between the 



true channel and the estimated channel. The channel estimation 
error variances can be written as 



= (i-aH 



(3a) 
(3b) 



B. Delayed feedback model 

There are two main selection strategies in the literature. In 
the first scenario, the relay nodes are in charge of monitoring 
of their individual source and destination instantaneous CSI 
0. In this strategy, CSI is estimated at the relay nodes for 
further decisions on which relay is the best node to cooperate. 
In particular, when the relay nodes overhear a ready-to-send 
(RTS) packet form the source terminal, they estimate the 
source link CSI, i.e., h em . Then, upon receiving a clear-to- 
send (CTS) packet from the destination, the corresponding 
destination link at the relays, i.e., h md is estimated. Right 
after receiving the CTS packet, the relay nodes ignite a timer 
which is a function of the instantaneous relay-destination CSfl 
The timer of the best relay expires first, and a flag packet is 
sent to other relays, informing them to stop their timers as 
the best relay is selected. Then the relays that receive the 
flag packet, keep silent and the selected relay participates in 
communication (6|. 

In the second strategy, the relay nodes estimate their individ- 
ual uplink CSI. The relays that correctly decode their received 
symbols from the source, send a flag packet to the destination, 
announcing that they are ready to participate in cooperation. 
The destination terminal, on the other hand, estimates the 
downlink CSI, orders the received SNR from each relay in 
the decoding set and feedbacks the index of the best relay 
that introduces the maximum received SNR via a log M bit 
feedback link. The selected relay then operates with full power 



In both strategies, a very important issue that must be taken 
into consideration is the selection speed. Communication links 
among terminals are time varying with a macroscopic rate in 
the order of Doppler shift, which is inversely proportional to 
the channel coherence time [6|. Any relay selection scheme 
must be performed no slower than the channel coherence time, 
otherwise selection is performed based on the old CSI, while 
the channel conditions are altered at the time that selection is 

'The timer is a function of relay-destination CSI for DF mode and source- 
relay-destination CSI for AF mode cooperation. 



performed. This might lead to a wrong selection of relay and 
therefore affect the performance of the system. 

In this paper, to simplify the notation, we adopt the second 
strategy in relay selection. However, the first strategy obeys 
the same rules and formulations. 

In the second strategy, since the feedback link only transmits 
the index of the selected relay, a lower feedback bandwidth is 
required. 

Let the estimated channel be h and the old estimated CSI 
based on which the selection is performed be h a . Since h and 
h are both zero mean and jointly Gaussian they can be related 
as follows [151 



^h + Jl-p f 2 v 



(4) 



where v ~ A/"(0, 1). We stress that h is the channel which 
is used for relay selection, whereas h is the channel which is 
used for decoding. 

III. Outage Probability of DF relay selection 

In the first time slot, the source terminal communicates 
with the relay terminals by broadcasting the signal x. Let the 
received signal at each relay be y sm , m = 1, . . . M, To decode 
the received symbol, each relay estimates the corresponding 
S —> R m channel and then decodes the received signal using 
an ML decoder. The old received signal at the m th relay can 
be written as 




message component error component noise component 

where we have substituted (l2bl into (f5ab . Using ( Bbl . the 
received effective SNR at each relay is given by 

PPe 2 \K m J 2 



where 7 

1 31 



7, 

fh 2 \h 



eff 



1 + Pol 



Pi 

1st 



(6) 



,l 2 



I — jr^i — . Similarly, in the second time slot, the 
effective SNR received by the destination via the m relay is 



imdp 



1 + Pol 



Pl^ 



(7) 



where 7 m 



ft 2 |/i | 2 

= l+ p^2 — , given that m € T>(s). Here, we as- 

U md 

sume that 7 s and 7 mdp are exponentially distributed with pa- 
rameters X sn J= (l + PalJ/p e and = {1 + Pa^J / p e . 

If the relay R m is in the decoding set, by sending a flag 
packet, it signals its capability of participating in cooperation. 
Then, based on the old channel realizations, the destination 
selects the relay with the best R m — > D link 



R 



,„. = arg max {7 }. 

m£T>(s) 



(8) 
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A. Outage probability with feedback delay 

Once the best relay is selected by destination, the index of 
R m * is fed back to all relays via a delayed feedback link. 
This means, at the time the relays receive the index, the 
system's CSI has changed, due to the time varying nature 
of communication links. Let the current R m * — > D channel 
realization, the corresponding estimate and the received SNR 
from the selected relay be /i m , d , h mM and j m * d , respectively. 
Due to feedback delay, relay selection is done based on 
past channel coefficients instead of current coefficients, i.e., 

7 m *d = 7md Where 171 = ar § maX {%«,}• 

i e T>(s) 

Assuming that the communication between source and 
destination targets an end-to-end data rate R, the system is in 
outage if the S — > R m * — > D link observes an instantaneous 
capacity per bandwidth^ C* = | log (1 + P"J mi , d ) that is below 
the required rate R, i.e., 



Pn 



Pr 



log(l + P%J<R 



x Pr 



l m * d < Po,7 m ^ = max {%^ }\T>{s) 



(9) 



x J2 Pr 

me~D(s) 



,(10) 



where R Q = — p — . 

B. Probability of decoding set T>(s) 

The relay, R m , is in the decoding set T>(s) if the S — > R Tl 
link observes an instantaneous capacity per bandwidth C sm 
that is above the required rate R 



C smp = -log 2 {l + P\ mp )>R. 



(ID 



Noting that \ is exponentially distributed, relay R m is in 
the decoding set if ifTOl 

Pr [R m G V{s)] = Pr[\ mp >R }. 

= exp (-\ mp R ) . (12) 

Finally, the probability of selecting a specific decoding set is 

m ' 

Pr[V(s)} = [] exp (-\ mo R ) 

meTi(s) 

x n [1 - exp (~\ smp R )] . (13) 

m<fD(s) 

2 Since the system is exposed to channel estimation error, the instantaneous 
capacity which is mentioned here, is only a lower bound for the true 
instantaneous capacity; hence, the derived expression for outage probability 
is a lower bound for the true probability of outage in the presence of channel 
estimation error and is exact with the perfect CSI. For further details please 
refer to Sec. |Vl] 



C. Outage probability conditioned on the decoding set T>(s) 

Conditioned on the decoding set with the old channel 
realizations, the outage probability with the new CSI is 



pV{s) = p r 

Defining 



l md < Ro,l md/J >max{7 i(lo }|:D(s) 



(14) 



X„ 



max {l idf }, 



e»(») 



and noting that \ m an d are independent, the conditioned 
outage probability on the decoding set is 

/>oo 

Po [s) = / Pr[l md <Ro\D{s)n mdp ] 
Jo 

x Pr [ 7mtto > Xm |2?(«), l mdfi ] f, m p m J d %n±o 
F (Ro)F(*, m Jf. d, ^ . (15) 



Conditioned on j m . and using |@), j md has a non-central Chi 
square distribution with two degree of freedom and parameter 

V m = c m 7 m dp' where 

write F (x) as ||26ll 



f h „ 



-b. Therefore we can 



k=0 



7(fc + l,^) 



(fc!) 



H2 



(16) 



where is the incomplete gamma function defined as 

j{s,x) = I r 

Jo 

Furthermore, F (x) is given by 

^» = n I 1 



(-A.j x) 



it > 



(_1)P>(«)I e -(^+\- 4 o+- A !£to ) :E . (17) 



Using (fTBI l and ( fTTb in ( fT3T > and noting that (a;) 



A e ^d." 31 and also the fact that 

mtLp 



x k exp(— ax)dx = 



k\ 

,fc+i ' 



the outage probability conditioned on the decoding set is given 
by ([T8T > at the top of the next page. 

Finally, substituting ( [13] ) and (ITSb into ( flOb . we obtain ( fl9l ) 
given at the top of the next page. 

For the special case where all the channel estimation errors 
and delayed feedbacks are the same, i.e., A = A = A^ = 
A idp , and c\ = c% = . . . c\x>( s )\ = c , (fl9l ) reduces to d20b given 
at the top of next page. 
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p v (s) y \l +1 7(fe+ l,A,„ d fl Q /2a 2 ) , c sk 
(fc!) 2 \ 2 J I 



■•+(-1) 



,ief(.) 



i. j, • ■ ■ .1 EB(') 

,1 



„ w ^ 7(fc + l,A md fl /2a 2 ) /ex* 
^ (fc!) V 2 J 



A^ 1 v 



E 



^ (A . + A.„ + A.. + A . 



-(-i) |cwl E 



^ (A d +A. d +%-) fe +! 
A..^ 1 



i,i i ex>( 

i # j , i 



(A m£to + A !(lo + A^ + . . . A ;do + A^)*+i 



(18) 



P °= E II exp(-A_i? ) [] [l-exp(-A_i? )] £ £ 



7 (fc+l ! A md i? /2 ( j 2 ) /c \* 



(fc!) 



ft) 



\ fc+i 



A^d, 



fe+1 



£™.U'+i + E 



+ A Ill£lo -f-) A;+1 . eD(ii) (A m£ i + A !ip + -^) k + 1 ^ ^ eI , w ( Am*. + A !tto + A 3dp + A^-f 4 -) 



.(_1)|T>WI ^ 



AJ 



fc+i 



,j, •■• ,l€I>( 

* i ( 



^Kndfi + A ido + A jtto + . . . \ ldp + A^-f-) 



fe+1 



(19) 



M 



P = V exp (-l\R ) [1 - exp (-Ai? )] M ' x i 



E 

fc=0 



7(fc + l,Ai? /2a 2 ) /c\* 



(fc!) 



(I) 



/ - 1 

TO — 1 



c\fc+l ' 



(20) 



IV. Average Symbol Error Rate 



In this section we derive a closed form expression for the 
ASER of the DF selection scheme considered in this paper, in 
the presence of channel estimation error and feedback delay. 
Let the received SNR via the selected path denoted by 7^ = 
Pj m t d - Then, the ASER can be written as 



Pe = aQ (^W\7 d ) f^ d (l m Jd^ d (21) 

where a and /3 depend on the modulation scheme, and fa t is 
the probability density function (PDF) of 7m „ d , which is given 
by (see Appendix for details) 



M 



+ E 



T>(s) 
|U( S )| > 1 



n a - n Bi 

itT>(s) i£V(s) 



(22) 



where S(x) is the delta function 



and Bj 



2 



1 



0P 



/3P + 2A fl 



(23) 



Here, f~ ._. .(a;) is the conditional PDF of the received 
signal via the selected path, which is given by 



(24) 



Noting that 

by inserting ( fT8l into (l25T l and the result in d24b the conditional 
PDF of the received signal is then given by (l26l i at the top of 
the next page. Finding a closed form formula for the integral 
in Pil l is not tractable. However, by substituting (l22l) into J2TT> . 
we obtain WIS at the top of the next page, where we have used 
the approximation | ]27l 



Q(x) 



E< 

n=l 



where 



(-1)™ +1 (A)™ 



(28) 



(29) 
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fe=0 



2a 2 k\ 



x fc+i 



Aj^J -1 



+E- 



fc+1 



(A . + A.„ + X.. + X 



. ( _i)P>wi 53 



Al +1 

771(2,0 



^ (A , +A +A , + ...A + A , 



(26) 



i=l u(») 

|U(s)| > 1 



Ha so n s « 



x a 



E EE 



» a n r ( & 



(n+1) 



^ V 2 / I 2<r 2 P/3 / 



)k=0 n=l I 1 ~r 2 p ( 3o-2 J 



(fc+1) 



\ fe+1 

An* \ "* 



(AjT 1C io A^^ 2 ) 



A fe+1 

+V- 



A^ 1 



+( _1)PWI 53 



5* J 



Pw + A^ + A,. + A. 



A^ 1 



,3, ■ ■ • ,1 £»(>) 
7t 3, I 



vrAhndp + ^ido + A Jdi0 + • • • A ;do + A md _ o -™-) 



fe+1 



(27) 



A/ 



p e = X+E 



z=i 



M 



(l-B)'S 



l uM-l 



x a 



; -11 „-l (1 + 2P^(T 2 ) fe+ 



I - 1 

x V m — 1. 
< V (30) 



i (™ + f ) 



£\t+l 



with ,4 = 1.98 and £ = 1.135. 

In the special case where all the channel estimation errors 
and delayed feedbacks are the same, i.e., A = \ = A mdo = 
X idf> , and Q = C2 = • ■ - Cmu)! — c > *E3 reduces to d30b at 
the top of the page, where B = B\ = . . . B 



V. Asymptotic Diversity Order 



decoding set, as 



PeWis) = « E E '" V A llH+il (f ) 
fe=0n=l I 1 T 2P/9CT 2 I 

l-l 

A \ x-^ \m — 1 



E 



2d 2 P/3/ ^ ( m +£)/=+!' 

771=1 v ^ 7 

g(n,k) ( X 

fe=0n=l l 1 + 2Po- i p> 



^EE 



fc+i 



(31) 



In the previous section, we have derived an exact ASER 
expression with feedback delay and channel estimation errors, 
which is valid for the entire SNR range. To gain further 
insights into the system's performance, we focus here on 
the high SNR regime and analyze the asymptotic ASER for 
the selection scheme under consideration conditioned on the 
decoding set, i.e., P e \V (s0 We assume that all the channel 
estimation errors and feedback delay parameters are the same, 
without loss of generality. Therefore, we have A = \ = 
A»*, = A ldo , and c x = c 2 = . . . C|i>( a )| = c 

Using d26b . we can write the ASER, conditioned on the 



3 In the DF scheme, selection relaying exploit the full diversity order, which 
is equal to the number of relays regardless of the cardinality of the decoding 
set |8|— 1 10]. Hence, in this work, we focus on the analysis of the asymptotic 
ASER conditioned on the decoding set 



where g(n, k) is a constant depending on k and n. 
In the following, we consider two different scenarios: 



A. Diversity order with feedback delay and perfect CSI 

With perfect CSI, we have A = 1. Therefore, as P — > oo, 
the dominant term in OTb would be the term corresponding 
to k — and n = 1, i.e. 



1 



e\V(s) 



1 + k /P 
1 



1 

~P' 



P + k 

oiP- 1 ) 



(32) 



Thus, in this case, the achievable diversity order in presence 
of delayed feedback is 1. 
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B. Diversity order in the presence of feedback delay and 
imperfect CSI 

In the presence of channel estimation error, since A is 
proportional to P, it should be taken into consideration. From 
OTb . we have 



lim P e ip( s ) 



OC 



lim 

p— »oc 



EE 



,k=0 n= 



= K. 



(33) 



Hence, the asymptotic diversity order in the presence of 
feedback delay and channel estimation errors is 0. This means 
that the outage probability or the ASER curves are expected 
to be saturated at high SNR with imperfect CSI. 

VI. Average capacity 

In this section, we derive an analytical expression for 
average capacity in the presence of channel estimation error 
and feedback delay. At the destination and after matched 
filtering, the received signal over the selected link is given 
as 



Ph* m „ 
Ph* 



= VPh 

= PPe\K. 



Ph* 



- u t 



(34) 



message component error component noise component 



Defining 

q = Ppk\K^\^ 

: 1 

(l34l can be written as 



n = Ph* 



(35a) 
(35b) 



x = qx + n. (36) 
The capacity conditioned on the decoding set is defined as 

C v ^ = arg max J(x,x) 

= arg max {H(x) — H(x\x)} 
= arg max {H(x) — H(n\x)} , 



where f x (x) is the PDF of the input signal x and %(■) is the 
differential entropy function. It can be deduced easily from 
d36i > that h and x are uncorrected, i.e., Ejria;*} = 0. However, 
n and x are not independent and therefore H(n\x) ^ H(n). 
This means that the standard procedures for deriving capacity 
expressions are not applicable to d37b . Alternatively, noting 
that conditioning does not increase the entropy, i.e., 

U{n\x) <H(n), 

we can lower bound the capacity as 



C v{s) > arg max {H{x) - H{h)}. 



(38) 



It must be noted that the distribution of n is not Gaussian. 
However, following PP . Theorem 1, we may assume that h 



has zero-mean complex Gaussian distribution with the same 
variance, which is the worse case distribution. In this case, we 
can re-write (f38T> as 



C v{s) > arg max {U(x) - U(n)\ > arg max { 

/*(*) U(x) 

H{x)-H{rJ)}, (39) 

where nt - CN (o, P 2 \h m „ d \ 2 (T 2 mtd + P\h mM | 2 ) . Maximizing 
(f39b with respect to f x (x) and using d36l ), a lower bound on 
the average capacity, conditioned on the decoding set, can be 
written as 

C V{ S) > C V{S) 



(40) 




= -log(l + P7 m *J 



Pe 2 I h I 2 

where j = = . D m o d — . Hence, a lower bound on the S-DF 

<m*d l-\-Pa^ ' 

capacity can be obtained as 

a h = J2 c ^ s) Prp(s)}- (4D 

Therefore, the average capacity bound, i.e., C\h, can be written 

as 



T>(s 



x / log (1 + P\ M ) U (%, d )d^ d . (42) 



Substituting d26] > into (l42l) yields d44l given at the top of 
the next page. Note that in our derivation we have used the 
integration formula [30] 

k\ 



/ ln(l + ax)x k e~ x dx = 2J jr — 

J n V 



{ > . ml-- 

ar \ a 



k—(J, 



! -- 



k — fJ, — t 



(43) 



(37) where 



/>oo y 

Ei(7) = / — exp(— x)dx 

J-y x 



is the exponential integral function. 

In the special case where all the channel errors and feedback 
delay parameters are the same, the average capacity bound can 
be obtained as in (l45T l at the top of next page. 

VII. Simulation Results 

In the section, we investigate the performance of selection 
cooperation in the presence of channel estimation errors and 
feedback delay through Monte-Carlo simulation. The transmit- 
ted symbols are drawn from an antipodal BPSK constellation, 
which means, a = 1 and (3 = 2. The node-to-node channels 
are assumed to be zero mean independent Gaussian processes, 
with variance <r? = a\ = 1 — a 2 e . In this case, h and 
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C] 



lb : 



In 2 



E 



II exp(-A_i? ) J] [l-cxp(-A sm , oJ R )] 



oo k 



md 

2^P 



k—fi 



A 



2a 2 P 



k—fi—t 



E\ " \ " ( C 2*) 
(k- u)\ 



\ fc+1 

£™.U'+i 



- - : eX 1 ^' 



kS 1 



(A^ + A^) 



(A mdo + A !4o + A mdo ^+i 



\ fc+ 



fc+i 



+ E (a — A. + A - - A , -r-i 



md,o 



^ (A , + A. , + A . + . . . A, , + A , 



fc+1 



(44) 



M 



Cib = E exp ( _ ^i? ) [1 - exp (-Ai? o ; 



oo k ( C ) 

I ' X|X EE A*+i(fe-//)! 



(- i )"r i e(^) E? 



k — f^ — t 



E 



/ - 1 

m — 1 



- (m + |) 



T - (45) 



/i are zero mean Gaussian processes with unit variance and 
p e = 1 — o"g. The variance of noise components is set to 
N = 1 and i? Q = 1 bps/Hz. 

Outage Probability: Fig. [2] shows the performance of the 
system with M = 4 relays with feedback delay (FD) for p. = 
0.6,0.7,0.8,0.9 and 1. In Fig. [2] we assume perfect CSI, i.e., 
p e = 1- We observe a perfect match between the analytical 
and simulation results. It can also be deduced from the slope 
of the curves that, for M = 4 and p f < 1, the diversity order 
of the selection scheme for is 1 . In the case of ideal feedback 
link, i.e., p f = 1, the diversity order of 4 is observed. 

Fig. [3] illustrates the performance of the system for M — 
2,3,4, with p. = 0.9,1. In Fig. [3] the same slope for 
different number of relays is noticed, confirming our analytical 
observations. For the case of p , = 1 and choosing M = 3,4, 
the diversity orders of 3 and 4 are observed, respectively. 

In Fig. |U we study the effect of both feedback delay 
and channel estimation errors. An error floor, in presence of 
channel estimation error, is noticed, as predicted by (36). 

ASER performance: Fig. [5] shows the ASER in the presence 
of feedback delay for M = 3. We assume perfect CSI 
knowledge. Assuming an ideal feedback link, the full diversity 
order of 3 is achieved. However, with feedback delay, the 
diversity order is reduced to 1, as predicted earlier. 

Fig. [6] illustrates the performance of the ASER in presence 
of channel estimation error and feedback delay for M — 2. It 
is clear from Fig. [6] that channel estimation errors reduce the 
diversity order of the system to zero, confirming our earlier 
analysis. 

Diversity order: Noting that the asymptotical diversity order 
d is given by the magnitude of the slope of ASER against 
average SNR in a log-log scale J2): 



d= lim 

SNR->oo log SNR 



(46) 



in this subsection we analyze the asymptotical diversity order 
of the underlying selection scheme. We assume perfect CSI 



knowledge, unless otherwise indicated. 

Fig. [7] shows the diversity performance of the system in 
presence of feedback delay link for different values of p f . 
We assume perfect CSI knowledge. It is observed that the 
asymptotical diversity order of the system tends to 1 for all pf 
values. This illustrates the destructive effect of feedback delay 
and demonstrates that relay selection in this case is annihilated. 

Fig. |8] illustrates the asymptotical diversity for different 
number of relays i.e., M = 2,3,4. It is obvious that at high 
SNR the diversity order is independent of the number of relays 
in the system. 

Fig. [9] depicts the asymptotical diversity order in presence 
of channel estimation errors. It is noticed that the asymptotical 
diversity order in this case is reduced to zero, confirming our 
earlier observations in Figs. [4] and [6] 

Average Capacity: Fig. \W\ shows the lower bound average 
capacity in bits per second per Hz per bandwidth versus 
SNR. It is observed that the presence of channel estimation 
errors result in capacity ceilings in the average capacity 
curves. It can be also seen that feedback delay aggravates the 
average capacity performance of the system. However, channel 
estimation errors have a greater effect in worsening the average 
capacity performance of the system. 



XI. Conclusion 

In this paper, we discuss a relay selection scheme in DF 
networks. We show that the presence of channel estimation 
errors and feedback delay degrades the performance and also 
reduces the diversity order of S-DF. We derive an exact ana- 
lytical expressions for the outage probability, average symbol 
error rate and average capacity bound. 

Appendix 
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Fig. 2. Outage Probability for M = 4 and perfect CSI in presence of delay 
in the feedback link. In this figure P e = 1 and p f = 0.6, 0.7, ... 1. 




SNR(dB) 

Fig. 3. Outage Probability for M = 2, 3,4 and perfect CSI in presence of 
delay in the feedback link. In this figure p e = 1 and p f = 0.9. 




5 10 15 20 25 30 

SNR(dB) 



Fig. 5.ASER for M = 3 and perfect CSI in presence of delay in the feedback 
link. In this figure p c = 1 and pf = 0.6, 0.7, . . . , 1. 
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S-DF, with DF and CE error, simulation, M=2 

o S-DF, with DF and CE error, analytical, M=2 
S-DF, perfect CSI and ideal feedback link, M=2 




FD and CE Error 
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pe = 0.9 








e 9 1 

O O j 








p E = 0.95 

\p. = i : 














u 




ideal feedback link and perfect CSI 















5 10 15 20 25 30 

SNR(dB) 



Fig. 6. ASER for M = 2 presence of delay in the feedback link and channel 
estimation error. In this figure p { = 0.9. 
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SNR(dB) 



Fig. 4. Outage probability for M = 2 presence of delay in the feedback link Fig. 7. Effective diversity order versus SNR. p c = 1 and M = 4, The figure 
and channel estimation error. In this figure p^ = 0.9. is plotted for different values of pf and is obtained from (46). 
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Deriving the Probability density function of 7^ 

Let 7 m<d denote the normalized received SNR at the desti- 
nation terminal over the S —> — >• D link. Then, the PDF 
of 7 m « d is written as 

f-v (x) = Prfall relays are off) f I , (x) 
J7 m *d V ; V ' ;J 7 mt<J | (all relays are off) V ; 



Fig. 8. Effective diversity order versus SNR. p e = 1, p^ 
2, 3, 4. The figure is obtained from (46). 
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— S-DF, average capacity, simulation, M=2 




o S-DF, average capacity, analytical, M=2 




A- = 1 




p/ = 0.7,0.8,0.9,1 




PS = 0.9, 1 




X^^fS Pe = 0.85,0.9,0.95 














f* Dashed line: average capacity with ideal feedback link 
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SNR(dB) 
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Fig. 10. Lower bound average capacity versus SNR for different values of pf 
and p e - 



^ ' Pr (relays in T>(s) are on) / 



|f(»)l > 1 



%n*fi I ( re ^ a y s i Q ^(s)are on) 



(x). (47) 



The i* 71 relay decodes its received signal erroneously with 
probability Bi which is given in (f23t . Since all S — > i?i links 
are statistically, we have 



50 100 150 200 250 300 350 400 450 500 
SNR(dB) 



Pr (all relays are off) = Bi. 



(48) 



On the other hand, if all the relays are off, then no communi- 
cation would occur between source and destination terminal. 
The received SNR at the destination terminal would be zero. 
Therefore, the conditional PDF can be written as 



Xx) = 5(x). 



(49) 



The probability of decoding set, given that there is at least 
one relay in T>(s), is given by 



Pr(relays in D(s) are on) 



50 100 150 200 250 300 350 400 
SNR(dB) 



Fig. 9. Effective diversity order versus SNR. p^ = 0.9 and M = 3, The 
figure is plotted for different values of p e and is obtained from d46t . 



Y[(l -Bi) [] B i ■ <50> 

ifV(s) ieV(s) 

Inserting (l48l . (|49l and d50]l in $7$ yields (l22l l. 
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